Generalized thermo vacuum state derived by the partial trace method 
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By virtue of the technique of integration within an ordered product (IWOP) of operators we 
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I. INTRODUCTION 

In nature every system is immersed in an environment, 
the problem about the system interacting with the en- 
vironment is a hot topic in quantum information and 
quantum optics. To describe system-environment evo- 
lution Takahashi-Umezawa introduced thermo field dy- 
namics (TFD) [l], with which one may convert the 
calculations of ensemble averages at finite temperature 



{A)=tr{Ap)/Z{P), p^e 



(1) 



to equivalent expectation values with a pure state |0(/3)), 



I.e., 



{A) = {0{l3)\A\0m, 



(2) 



where P — l/kT, k is the Boltzmann constant, and 
Z {(3) = trp = tre^^^ is the partition function; H is 
the system's Hamiltonian. Then how to find the explicit 
expression of |0(/?))? If one expands |0(/3)) in terms of 
the energy eigenvector set of H, |0(/3)) = ^„ \n) fn{(3), 
and then substituting it into Eq.([2]), which results in 
fl{(i)fm{li) = Z-^ (/3) e-f^^-5n,m (after comparing with 
Eq.©). By introducing a fictitious mode, (n |m) = 5n,m, 
then Takahashi-Umezawa obtained 



|0(/3)) = 



(3) 



which is named thermo vacuum state, and S {9) thermo 
operator. 



S (9) EE exp [9 



a^a^ — aa)] , 



(5) 



which is similar in form to the a two-mode squeezing 
operator except for the tilde mode. is a parameter 
related to the temperature by tanh^ = exp {—^p) ■ 
An interesting question thus challenges us: For the 



t2 



then what 



Hamiltonian being H = uia^a + K*a 
is the corresponding thermo vacuum state? One may 
wonder if this question is worth of paying attention since 
this H can be diagonalized by the Bogoliubov transfor- 
mation as a new harmonic oscillator, correspondingly, the 
thermo vacuum state for H can be obtained by acting the 
same transformation on |0(/?)) in ^ (see Eq. (A9) in the 
Appendix). To make this issue worthwhile, we empha- 
size that we shall adopt a completely new approach to 
construct thermo vacuum state and our result is simpler 
in form than that in Eq. (AlO). Our work is arranged as 
follows. In Sec. 2 by re-analyzing Eqs. ([I])-© we shall 
introduce a new method (the partial trace method) to 
find the explicit expression of |0(/3)) in Q. Then using 
this method, we obtain the expression of |0(/3)) in Eq. Q 
in Sec. 3. For the degenerate parametric amplifier, we 
derive a generalized thermal vacuum state \4>{f3)) in Sec. 
4. Section 5 is devoting to presenting some applications 
of |0(/3)). 



Thus the worthwhile convenience in Eq.® is at the ex- 
pense of introducing a fictitious field (or called a tilde- 
conjugate field, denoted as operator a^) in the extended 
Hilbert space, i.e., the original optical field state \n) in 
the Hilbert space Ti. is accompanied by a tilde state \n) in 
Ti.. A similar rule holds for operators: every Bose anni- 
hilation operator a acting on Ti. has an image a acting on 
Ti., [a,a^] — 1. These operators in Ti are commutative 
with those in Ti. 

For a harmonic oscillator the Hamiltonian is huja^a, 



/Vnl |0) , Takahashi-Umezawa obtained the ex- 
plicit expression of |0(/3)) in doubled Fock space, 

|0(/3)) = sech6'exp [a^a^ tanh 6l] |o6) ^ S (9) |00) , (4) 



II. THE PARTIAL TRACE METHOD 

Following the spirit of TFD, for a density operator p = 
e^^^ / Z [P) with Hamiltonian H, we can suppose that 
the ensemble averages of a system operator A may be 
calculated as 



A^tr{pA) = {ij{(3)\A\^{(3)), 



(6) 



where |V'(/3)) corresponds to the pure state in the ex- 
tended Hilbert space. 

Let Tr denote the trace operation over both the system 
freedom (expressed by tr) and the environment freedom 



2 



by tr, i.e., Tr = trtr, then we have 

(^(/3)|A|V(/3)) - Tr[A|V-(/3))(V(/3)|] 

= tr[yltr|^(/3))(V.(/3)|] 



Note that 



(7) 



(8) 



since involves both real mode and fictitious mode. 

Comparing Eq.(I7]) with Eq.([T]) we see 



tr|^(/3)) (V'(/3)| 



IZ{(3). 



(9) 



Eq.® indicates that, for a given Hamiltonian if we 
can find a density operator of pure state in dou- 

bled Hilbert space, whose partial trace over the tilde free- 
dom may lead to density operator e~^^ /Z {(i) of the sys- 
tem, then the average value of operator A can be calcu- 
lated as an equivalent expectation value with a pure state 
|V(/3)) , i.e., {A) = tr {Ae-^" /Z (/?)) = {^{f})] A 

In particular, when H = huia^a, a free Bose system, 
Eq. ([9|) becomes 

ti |0(/3)) (0(/3)| = (1 - e-P'''^) e-^^^"'''' = p„ (10) 

Pc is the density operator of chaotic field. This equa- 
tion enlightens us to have a new approach for deriving 
|0(/3)) : |0(/3)) (0(/3)| in doubled Hilbert space should be 
such constructed that its partial trace over the tilde free- 
dom may lead to density operator pc of the system. In 
the following we shall employ the technique of integration 
within an ordered product (IWOP) of operators [3, 0, [E] 
to realize this goal. 



III. DERIVATION OF 0(/3)) IN EQ.® VIA THE 
NEW APPROACH 

Using the normally ordered expansion formula [6j 



exp { (. 



(11) 



(where the symbol : : denotes the normal ordering form 
of operator), and the IWOP technique we have 



: exp{(e-'3'"'"-l)ata}: 



■ : e 



+zae 



(12) 



Remembering the ordering form of vacuum projector op- 
erator |0) (0| =: e""^": , we can rewrite Eg. p^ as 

: exp{(e-'3'^"^- l)a^a}: 
where | z) is the coherent state 0, 01 in fictitous mode 



\z) ~ exp i^zd) — z*a) |0) ,a\z) ~ z \z) , (6 \z) = 



(14) 



Further, multipling the factor (l — e ^^'^) to both sides 
of Eq. (fT3|) and using the completeness of coherent state 
/ ^ 1^) {z\ = 1, we have 



1-e 



(1 



X Eq.dl 



1-e 



a z 

TT 
TT 



(z| e^*"'^"""^'' |00)(06 



00) (oo| 



= tr[0(/3)(0(/3)|] 
where 



5) 
(15) 



|0(/5)) = VT^TT^ 



exp 



|00), (16) 



which is the same as Eq.(|3]). Thus, according to Eq. pH)) . 
from the chaotic field operator we have derived the 
thermo vacuum state, this is a new approach, which has 
been overlooked in the literature before. 



IV. GENERALIZED THERMO VACUUM STATE 

Now, we consider a degenerate parametric amplifier 
whose Hamiltonian is 

H = uja} a + K* a)'^ + KO^ , (17) 

whose normalized density operator p is 

p (tre-'^^) = e-P" = ^-f^^-^ , (ig) 

Recalling that ^ (a^a -|- i) , ^a'^^ and ^a? obey the 
SU(1,1) Lie algebra, thus we can derive a generalized 
identity of operator [l3] as follows: 

exp [/a^a + ga)"^ + /ca^] 

where we have set T>^ — f ^ — '^kg. Thus Comparing 
Eq.([TH]) with Eq.lll]) we can recast Eq.dlHl) into the fol- 
lowing form 



(tre"''^) p = exp [E*a'''^] exp [a^a In A] exp [Ea^] , 

(20) 

where we have set 



D 



uj sinh I3D + D cosh I3D ' 



(21) 



E = -^KsinhpD. 
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Further, using the formula in Eqs. (fTT|) and ((T3l) . we have 
{tre-^") p 

= VAe^" exp [E*a^^] : exp {(A - 1) a^a} : exp [Ea^] 
= VX^ J ^6^*'^''+^^*''' |0)(0|e^'^'+^^"(z |0)(0 



^/A^ I — (z| eE'a^'-'+Vxz'a'' |qq^ ^qq giSa^+VAza | ~^ the idea from Eq.H]) to (H]), the system operator A can 

be calculated as (A)^ = (0 {P)\A |^ (/?)) . Thus uisng the 
completeness of coherent state and the integral formula 
Eq.dlSl) as wen as noticing (0 {f3) {(3)) = 1, (1 - A)^ - 



V. APPLICATIONS OF GENERALIZED 
THERMO VACUUM STATE 

A. Internal energy distribution of the system 

^) As an application of Eq. ((27)) . we can evaluate the each 
term's contribution to energy in Hamiltonian. Based on 



Ta^ 



tr 



1 00) (00 



-y/Xa 



tre 



)tr[|</)(/3))(0(/3)|], 



which indicates that the pure state in doubled Fock space 
for Hamiltonian in (|17p can be considered as 



l0(/5)> 



Z{f3) 



00) , (23) 



(22^ \Ef = 4Asinh^ {f3D/2) , then we have 

(ioa^a)^ = c.(0(/3)|(aat_i)|^(/?)) 
= 2LuX^/'^smh{l3D/2)-^ 



where the partition function Z {(3) is determined by 
Z iP) = tre-^" = tr { VA^e^*"^'' : e^^"!)'^''^ : e^"^' } . 

Using / - 



(24) 



z) {z\ = 1 and the integral formula [lll | 
exp (^(\zf+^z + r^z* + fz^ + gz*^^ 

whose convergent condition is Re(C ± / ± (?) 
< 0, we can get 



(25) 



< 



Re( C'~ifa] 



(l-A)'-4|£;r 2sinh(/3D/2)' 



(26) 



Thus the normalized state for Eq.([T8|) in doubled Fock 
space is given by 



10 (/3)) = y'2Ai/2sinh(/3i:)/2)e^ 
and the internal energy of system is 

which leads to the distribution of entropy 



|00), (27) 
(28) 



S = -ktT{p\np) = ^{H) 
D 



fclnZ(/3) 



= — coth {f3D/2) - fc In [2 sinh {PD /2)] . (29) 



In particular, when k = 0, leading to D — uj, so Eg. (1271) 
just reduces to |0(/3)) with uj huj, and Eas. ipS)) and 
(pg]) become | (coth(/3w/2) - 1) and ^ coth (/?tj/2) - 
fcln[2sinh(/3a;/2)] , respectively, as expected [l^- Thus 
by virtue of the technique of IWOP we can display the 
partial trace method to deduce the pure state representa- 
tion for some new density operators of light field at finite 
temperature. 



d^z 



-(l-\)\z\^+Ez^+E'z' 



Ai/2 



d 2tJsinh(/3L>/2) 



(l-A)^-4|i;|^ 



|(^coth/37^/2-l 



(30) 



and 



dE* 



{l-Xf 



= -■Ll-coth/3D/2, 



as well as 



(ko^)^ = -l±-cQth.pDl2. 



4EE* 
(31) 

(32) 



From Eqs. ((3Tl) and (|32|) we see that the two items {K*a^^ 
and Ko^) have the same energy contributions to the sys- 
tem, as expected. Combing Eg s. ([50)) - ([5^ we can also 
check Eq. (| 



B. Wigner function and quantum tomogram 

The Wigner function plays an important role in study- 
ing quantum optics and quantum statistics [l^, . It 
gives the most analogous description of quantum mechan- 
ics in the phase space to classical statistical mechanics 
of Hamilton systems and is also a useful measure for 
studying the nonclassical features of quantum states. In 
addition, the Wigner function can be reconstructed by 
measuring several quadratures P {xg = x cos 9 + p sin 9) 
with a homodyne detection and then applying an in- 
verse Radon transform — quantum homodyne tomogra- 
phy [l^. Using Eg. (|^5|) one can calculate conveniently 
the Wigner function and quantum tomogram. Recalling 
that the single-mode Wigner operator A (z) in coherent 
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state representation is given by [T^ [T3| 



A (a) 



-z\ e 



2{za* ~z* a) 



(33) 



thus the Wigner function is 
W{a) = (0(/3)|A(a)|0(/3)) 



(0(/3) \z) i-z |0(/3))e- 



-2{za* — js* ck) 



^ tanh (/?j:'/2) ^..p. [^|„|2_^(-^^2_^^.^.2-)j ta„h(/3I^^^J 

where we have noticed (1 + A)^ — 4|£'|^ — 

u^sL^tiD+DcoiliiD and used the integral formula 
(l25)l . In particular, when k = 0, Eq. p4|) reduces to 



W{a) 



tanh (/9i:'/2) 



exp|-2|a|^tanh (/5i:'/2)| 



(35) 

which is just the Wigner function of thermo vacuum state 
|0(/3)). 

On the other hand, we can derive the tomography 
(Radon transform of Wigner function) of the system by 
using Eq. (f26|) . Recalling that, for single-mode case, the 
Radon transform of the Wigner operator is just a pure 
state density operator [18], 



(36) 



5 [q - fq' - gp') A (a) dq' dp' - \q) ^ ^ {q\ 
where a = {q + ip)/\/2, and (/, g) are real. 



\q) 



C exp 



V2 ^ 



qa' 



,t2 



A 2 

1/4 



|0) , (37) 



and C = [n{p+g^)] ''%xp{-gV[2(/2 + ^2)]}, ^ ^ 

f—ig— a//^ + g^e"^"^. Eg. ([57)) is named as the interme- 
diate coordinate- momentum representation [T^. From 
Eq. ([5S)) and Eg. ([^5]) it then follows that the tomogram 
can be calculated as 



, we ob- 



In sum, by virtue of the technique of integration within 
an ordered product (IWOP) of operators we have pre- 
sented a new approach for deriving generalized thermo 
vacuum state which is simpler in form that the result by 
using the Umezawa-Takahashi approach, in this way the 
thermo field dynamics can be developed. 

Appendix: 

As a comparison of our new approach with the usual 
way of deriving thermo vacuum state in TFD theory, 
in this appendix, we shall derive the explicit expression 
of |(/>(/3)) by diagonalizing Hamiltonian (|17p . For this 
purpose, we introduce two unitary operators: one is a 
single mode squeezing operator, 



S — exp 



t2 



v a 
M~2 



exp 



1 



In- 



1 



exp 



/i 2 



(Al) 

where fi and v are squeezing parameters satisfying the 
unitary-modulate condition — — 1; and the other 

is a rotational operator, R = exp (^^a^a^ , which lead to 

the following transformations, 

SaS^ = lia — iya\ So' = ^a) — va, (A2) 

S^aS = /ia + va^ , S^a^ S = ^jlo' -t- va, (A3) 

and 

i?ai?^ = ae"^, Ra) = a^e^. (A4) 

Thus, under the unitary transform SR^ we have (setting 

K = \K\e"l') 

H' = SRHR^S^ = ojSa^aS^ + \k\ Sa^^S^ + |k| Sa^S"^ 
— (uJiM^ + lliv^ — 4 |k| /ii/) o) a + (lov — 2 |k| /i) v 

In order to diagonalize Eq.(A5), noticing [P' — v^ = 1 and 
making |k| (/i^ + i^^) — LOjiv = 0, whose solution is given 
by 



1 



M = 



n {q)f^g ^ {<P (/3) {q \<j> (/3)) =J— \f,s ill \<P W) I 

(38) 

Then submitting Eqs.dS?]) and ^ into Eq, 
tain 



2w' ' 2'^ ~ 2^ 
then Eq.(A5) becomes 

2 



2'" 



4k| 



(A6) 



(A7) 



T^iq)f. 



2sinh(/3D/2) 



f.9 



i.e., the diagonalization of Hamiltonian is completed. 

According to Eg. (|16p . the thermal vacuum state corre- 
sponding to density operator p' = e^^^ /tr {^e^^^ ^ = 



C-\IX + \GY /\ 

, ^ 1 - AReG-i 

X exp < 2q 



—f3uj'a a 



is given by 



\Af(X+\Gnx 



Re 



2E 
AH] 



atate-'5-'/2 



where we have used Eq.(l25l) and set G = 1 -I- 2e"'^£;. 
Eg. ([55)1 is the positive-definite tomogram, as expected. 
As far as we are concerned, this result has not been re- 
ported in the literature before. 



|0(/3)) = Vl - e-P"^' exp 
Thus the generalized thermal vacuum state is 



00) . (A8) 



exp 



atate-^"'/2 



|00) 



(A9) 



5 



Using the transformation in (A3), (A4) and noticing 
Eq.(Al) as well as - = 1 we can finally put 



lg-(/3^'+»0)/2^t5t 



Eq.(A9) into the following form 



2^ 



.t2 



2/x 



^t2 



1 00) , 



(AlO) 



Comparing Eq.(AlO) with Eq.(l27l), we see that Eq. (l27|) 
is simpler in form than that in Eq.(AlO). 
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